— 7r transitions in Josephson junctions with antiferromagnetic interlayers 
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We show that the dc Josephson current through superconductor-antiferromagnet-superconductor 
(S/AF/S) junctions manifests a remarkable atomic scale dependence on the interlayer thickness. At 
low temperatures the junction is either a 0- or 7r-junction depending on whether the AF interlayer 
consists of an even or odd number of atomic layers. This is associated with different symmetries of 
the AF interlayers in the two cases. In the junction with odd AF interlayers an additional n — 
transition can take place as a function of temperature. This originates from the interplay of spin-split 
Andreev bound states. Experimental implications of these theoretical findings are discussed. 

PACS numbers: 74.45.+C, 74.50.+r, 75.50.Ee 



The developing field of superconducting spintron- 
ics subsumes many fascinating physical phenomena 
with potential applications that may complement non- 
superconducting spintronic devices^]. In addition 
there is an increasing interest in the novel proper- 
ties of interfaces and junctions of superconductors and 
magnetic materials 0, 0- An important example is 
the — 7r transition in superconductor-ferromagnet- 
superconductor (S/F/S) junctions where, depending on 
the temperature and the width of the ferromagnetic in- 
terlayer, the ground state of the junction is characterized 
by an intrinsic phase difference of n between the two 
superconductors 0, 0, • This may be important for su- 
perconducting digital circuits, and has been proposed as 
a possible basis for quantum qubits^, 0, 0, M, E]- The 
interplay of magnetic and superconducting order leads 
to interesting mesoscopic physical phenomena also at in- 
terfaces between antiferromagnets and superconductors. 
Recent theoretical studies have shown that a characteris- 
tic spin dependent quasiparticle reflection at the AF sur- 
face, the so-called Q-reflection, combined with Andreev 
reflection on the superconducting side, leads to novel low- 
energy bound states near such interfaces ^l]- These An- 
dreev bound states have important consequences for the 
associated proximity effect and can, for example, be de- 
tected in tunneling spectroscop y as sub-gap peaks in the 
resulting local density of states[12j. 

In this Letter we investigate effects of the Q-reflection 
on the Josephson current in S/AF/S tunnel junctions in 
s-wave superconductors. An enhancement of the low- 
temperature critical current in such junctions was found 
in the limit of a sufficiently small ratio m/t, where t and 
m denote the hopping matrix element and the antifer- 
romagnetic order parameter, respectively |ll| . As shown 
below, the parameter m/t becomes sufficiently small in 
tunnel junctions only under the condition (m/t) <C \[~D, 
where D is the transparency coefficient. We will find that 
in the opposite case (m/t) > a/D, which is also realistic 
for fabrication, the Josephson current in S/AF/S junc- 
tions exhibits new interesting properties. At low temper- 



ature T, the current-phase relation reveals either a 0- or 
a 7r-junction state depending on whether the AF inter- 
layer consists of an even or odd number of atomic layers, 
respectively. This atomic-scale thickness dependence dif- 
fers from that of S/F/S junctions where the period of the 
alternation between 0- and 7r-states strongly depends on 
the exchange field which controls the proximity-induced 
damped spatial oscillations of the pairing amplitude in 
the ferromagnetic metal. In contrast, for the S/AF/S 
junctions the — 7r behavior is a true even/odd effect re- 
lated to the difference in symmetries of the odd and even 
AF interfaces and the corresponding interface S matri- 
ces, as well as the spectra of the Andreev bound states. 
For odd (110) AF interlayers, the supercurrent displays a 
prominent anomaly with increasing T revealing another 
7r — transition. Wc will show that this remarkable result 
is a consequence of the interplay of spin-split Andreev 
bound states contributing with opposite sign to the total 
supercurrent. Lastly, we discuss possible experimental 
consequences of the effects in question. Our theoretical 
analysis includes both the self-consistent numerical solu- 
tions of the Bogoliubov-de Gennes (BdG) equations and 
a quasiclassical analytical approach to the superconduct- 
ing leads, allowing us to describe the AF interface fully 
microscopically with atomic scale accuracy, and interpret 
the results physically. 

Model. The Hamiltonian is defined on a 2D square 
lattice with superconducting and magnetic order 
parameters, and lattice constant a = 1: 

h= - *E a l^ + E( A 4 £ a+ H - c -) 
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Here, c| CT creates an electron of spin a on the site i, fi 
is the chemical potential and hi a = c^c^ is the particle 
number operator. The associated BdG equations are 
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The diagonal blocks are given by K,^ — —tS^ — /xdy ± 
aniiSij, where a = +1/ — 1 for up/down spin and Sij and 
5(ij\ are the Kronecker delta symbols connecting on-site 
and nearest neighbor sites, respectively. The off-diagonal 
block T>ij describes s-wave pairing T>ij = —A{Sij. The 
net magnetization Mi = | [(^if) ~ (^ii)] an( i ^ ne P a i r ~ 
ing amplitude Fi = (ciic^) are related to m t and Aj by 
rrii = UiMi and Aj = —ViFi. The coupling constants 
Ui (Vi) are site dependent and non-zero on (off) the L 
atomic chains constituting the AF interlayer. This sta- 
bilizes the staggered AF order on the interlayer and the 
superconducting order outside this region. We choose the 
x (y) axis to run perpendicular (parallel) to the interface. 
The calculations for planar tunnel junctions with crystal 
periodicity along the interface are reduced to a ID prob- 
lem by Fourier transforming along y. This introduces a 
crystal- vector component k y as a parameter and the BdG 
equations have to be diagonalized for each k y . 

The dc Josephson current j rr 
boring sites r = (x, y) and r' 
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where < k y < and f(E) is the Fermi func- 

tion. For the (100) junction x' = x + 1 and the current 
takes the form 0, but without the cos(k y /y/2) factor and 
with — 7T < k y < 7r. Since the vector potential generated 
by the current is not taken into account, we discuss only 
planar junctions in the tunneling limit when the current 
j is well below the thermodynamic critical pair-breaking 
current. Below, we fix the phase of the superconducting 
order parameter at each end of the system, obtain the 
self-consistent solutions, and calculate the current based 
on Eq. (0) . As is well-known, the current is only conserved 
when the superconducting order parameter is calculated 
fully self-consistently^^. For more details on the nu- 
merical and analytical appr oaches used in this Letter, we 
refer the reader to Ref. [l2|. 

Results. We have studied the dc Josephson current in 
both (100) and (110) S/AF/S junctions. Below, how- 
ever, we focus on the (110) interfaces where the effects 
in question are more pronounced. Spins along (across) 
a (110) AF interlayer are identically aligned (alternate). 
Fig. nji,b show two representative T = current-phase 
relations for the S/AF/S tunnel junctions with different 
thicknesses L of the sandwiched (110) AF layer. The 
curves display striking 0- or 7r-junction sinusoidal behav- 
ior depending on whether the AF interface contains an 
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FIG. 1: (Color online) Current-phase relation at T = 0, U = 
8.0t (a,b) and temperature dependence of the critical current 
3c = li(f )| (c,d) for S/AF/S (110) junctions with L = 6 (a,c) 
and L = 5 (b,d). In (c,d) the graphs have been normalized to 
their values at T = 0. (e) Oscillations of the critical current 
as a function of the AF thickness L (U — 8.0i). (f) Current- 
phase relation at various temperatures along the U — 8.0t 
curve in (d). Parameters used in all figures: V = 2.0t, /i = 0. 
Currents are in units of et/M y . 



even (Fig. [TJi) or odd (Fig. 2Jd) number of chains, respec- 
tively. Fig. ^ shows the critical current as a function of 
thickness L. The temperature dependence of the critical 
current j c (T) is shown in Fig. ^,d for varying magnetic 
strength U. Clearly, the magnetism leads to remark- 
able anomalous T dependence of the critical current. For 
small values of U, where the influence of the Q-rcflection 
is more pronounced since the Andreev bound states have 
lower energy, j c {T) deviates from the standard behav- 
ior by varying monotonically for the even junctions and 
non-monotonically for odd L. In particular, we find that 
in addition to the low T alternating — ir transitions as 
a function of L, the odd junctions can exhibit another 
7r — transition as a function of T . This n — transition 
is clearly shown in Fig. where we plot the current- 
phase relation at various T along the U = 8.0t curve in 
Fig. ^i. Note that in Fig. \Q,d the conventional behavior 
near T c , j c ~ [1 — T/T c ], is recovered only in the limit of 
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FIG. 2: Eigenbands for the S/AF/S junctions with U = 8.0t 
and fj, = T = ip = 0.0 and L = 6 (a) and L = 5 (b). As seen, 
the even or odd number of magnetic layers lead to qualita- 
tively different bound state bands inside the gap. 

large U . The anomalous temperature dependence near T c 
appears to be associated with the magnetism leaking into 
the superconducting leads by the proximity effect 12]. 

In junctions between identical superconductors with a 
thin interlayer L <C £, s ~ ^ v f/\^\, the whole Josephson 
current is carried through the interface by the phase- 
dependent subgap Andreev bound states. The bound 
states at the two separate S/AF interfaces [ill IT^ . mix 
in the junction geometry resulting in qualitatively differ- 
ent bound state bands for the odd or even interlayers. 
This is seen in Fig. [21 where we show the spin-down 
eigen-spectrum for a even (a) and odd (b) (110) S/AF/S 
junction as a function of k y . In the case of odd L only 
the total spectrum, including both spin-down and spin- 
up states, is symmetric with respect to the Fermi sur- 
face. For larger values of U, the transparency of the 
junction decreases, the specular reflection becomes more 
pronounced and the bound states move towards the gap 
edge and become more extended. 

We have carried out analytical calculations of the sub- 
gap spectrum under the conditions A <C m,t, L <C £ s > 
whereas our numerical self-consistent studies are ap- 
plicable in a wider regime. The normal-state spin- 
dependent reflection and transmission amplitudes are the 
same on both sides of a (110) odd AF interface: r a = 
y/Rodd exp(icr0/2), d a — io\JD a dd exp (iaQ/2). Here the 
transparency coefficient is spin- independent D dd(k y ) = 
cosh -2 [2(L+ i)/3], sinh/3 = m / [At cos(k y / y/2)} and 
Rodd(ky) = 1 - D od d(ky). The difference 8 = 
Of — &i between phases of spin-up and spin-down 
reflection amplitudes takes the form sin0(fc y ) = 

[m/2tcos(fcy/ % /2)] |l+ [m/4tcos(fc y /V2)] 2 | . If 

m < At cos(k y / y/2) , satisfies the relation tt/2 < < tt, 
whereas for m > At cos(k y / y/2) one gets < O < tt/2. 
The narrow regions of k y near k y = ±tt/^/2 turn out 
not to be important for the effects in question and in the 
estimations one can put cos(fc 1) /v / 2) ~ 1. 



The key properties of the odd S/AF/S junctions are as- 
sociated with the fact that they contain symmetric mag- 
netic interlayers with respect to the two superconductors. 
In the (110) odd interlayer both outermost chains belong 
to the majority spin polarization, defined as "spin-up" , 
and the Andreev states are spin-split. The physical ori- 
gin of the 7T — transition at finite T is just the interplay 
of spin-split Andreev bound states. The general struc- 
ture of the S matrix containing the present reflection 
and transmission amplitudes, is similar to that found in 
Ref. E3 

for symmetric ferromagnetic interfaces. An im- 
portant distinction to the AF case, however, is the very 
different expression for the parameter 0. The presence of 
the low-energy Q-reflection results in quite large values 
> tt/2, which ensure the tt — transition with varying 
temperature in a wide range of the AF order parameter 
m < At. Applying the quasiclassical equations to the 
superconducting leads and taking into account the AF 
interface properties within the iS matrix approach, we 
obtain the following spectrum for the Andreev bound 
states: Ef{k y ) = cr|A|sgn[sin (f ± S)] cos (f ± S). 
Here, sin S = \/D dd cos(x/2) and x is the phase dif- 
ference between the left and right superconducting leads. 
As seen from the energy spectrum E^(k y ), the parameter 
mj At becomes negligible in tunnel junctions under quite 
strict conditions m/At <C \/D dd- This implies ra tt 
when four dispersive bound states E^(k y ) reduce to two 
doubly degenerate states, leading to the particular results 
of RefE] for the Josephson current. 

Below we consider the different conditions, y/ D dd <C 
m/At, 1 which can be easily satisfied in tunnel junctions. 
The expression for the Josephson current which follows 
from the bound state spectrum in the case \JD Q d d *C 
m/At, 1, describes the tt — transition with increasing 
temperature for 7r/2 < < tt, and differs only by the ad- 
ditional factor of 1/2 from Eq. (4) of Ref. Q (with a = 1 ) 
valid for clean symmetric S/F/S tunnel junctions. At 
T = the current is given by a sum over k y of the expres- 
sions j°f y {T = 0) = -(l/2)e|A|L> odci cos(0/2), whereas 

near T c we get j°f y {T ra T c ) = -e\ A\ 2 D odd cos 0/4T c . 
From this we see that the product j°f y (T = 0)j°f y (T ra 
T c ) is negative for 7r/2 < < tt and positive for < tt/2. 
Hence, if < tt/2, there is no tt — transition with in- 
creasing T and the 7r-state remains the ground state of 
the odd (110) S/AF/S junctions for all T < T c . 

For the even (110) junctions, the outermost chains of 
the AF interface have opposite spin polarizations and the 
Andreev states are spin degenerate. In this geometry the 
spin-dependent reflection amplitudes on the two sides of 
the interfaces differ r% t(7 — —r* 2a = \/R ev exp(ier0/2). 
The transmission amplitude is real and spin-independent 
d = y/D ev = 1/ cosh [2(L + l)/3] and the expressions for 
and (3 coincide with the odd case. Such structure 
of the S matrix is characteristic also for the three-layer 
FIF interface with antiparallel orientations of the two 
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ferromagnetic magnetizations |15| . The corresponding 
spin-degenerate spectrum of the Andreev states in even 
S/AF/S junctions, E ev (k y ) = ±A[D ev (k y ) cos 2 (x/2) + 
Rev(k y ) cos 2 (0/2)] 1 / 2 , has the same form as that found 
for the FIF three-layer with antiparallel magnetizations 
(see Eq.(4) in Ref. [l5j), but with a significantly different 
expression for O. The Josephson current at zero tem- 
perature is a sum over k y of the expressions (T — 
0) = e| A\D ev /2 cos(0/2), whereas near T c we obtain 
3™k y ( T ~ T c) = e\A\ 2 D ev /4T c . There is no - ir tran- 
sition with varying temperature, in agreement with our 
numerical results and with the S/FIF/S junctions with 
antiparallel magnetizations. 

The different symmetries of even versus odd AF inter- 
layers are also responsible for the low T — ir transitions 
as a function of L. Indeed, the 7r-state is the ground 
state of clean S/F/S junctions with a = 1, 9 > 7r/20, 
whereas the 0-state is the ground state of the S/FIF/S 
junctions with antiparallel magnetizations [l5j . The di- 
rect analogy between the odd (even) AF interfaces and 
the F interfaces at a = 1 (the FIF interfaces with antipar- 
allel magnetizations), results in the correct sequence of 
the transitions at T = 0. The even/odd — ir transitions 
can also be related, within the perturbative approach, 
to effects of localized spin states in interfaces and the 
anticommutation of fermions[ll HI HI HI EI El] . 

The critical current j c in odd (110) S/AF/S junctions 
at T — is reduced by the factor cos(8/2) compared to 
the standard junctions with the transparency D a dd- For 
m/f < 1 we have cos(0/2) <C 1 and the relative suppres- 
sion is significant. The presumed condition \JD dd 
to/42; gives the smallest value cos (0/2) ~ \fD^dd and a 
reduced current j° dd oc D Z J d 2 d . Contrary to the odd case, 
the critical current j% v in even junctions at T = is en- 
hanced by 1/ cos (0/2) compared to standard junctions 
with the transparency D ev . The maximal relative en- 
hancement is realized for cos(0/2) ~ \JD ev , where the 
applicabilities of the present results border to those of 
Ref. El Then ~ e|A|-^/D el ,, in qualitative agreement 
with the net critical current at T = found in Ref. fill 

Fabrications of ultrathin interfaces with atomic-scale 
control of the thickness over a macroscopic area, similar 
to the case of ultrathin films [22], |23j, would allow ob- 
servations of the even/odd effect for S/AF/S junctions. 
Alternatively one should average the Josephson current 
over interface imperfections. For example, assume that 
the procedure can be reduced to an averaging over thick- 
ness variations within a few layers. If L 3> 1, we get 
Dodd{L) w D ev (L ± 1). For m<Ct the critical currents 
are related as jf ~ (t/m) 2 \j° dd \ > \j° dd \ at T = 0, 
whereas near T — T c they are of the same order and 
sign. Since the current jf strongly dominates \j° dd \ at 
low T, the result of the averaging is that the ground state 
is the 0-state with the anomalous critical current j™ (T) 
and there is no — 7r transition, if to <C t . For m ~ t even 



and odd currents are of the same order and therefore sim- 
ilar samples can have differing signs of j c . If < 7r/2, 
jf and j° dd have opposite signs and we expect a pro- 
nounced suppression of the net critical current at all T. 
For > 7r/2, the odd and even currents have opposite 
signs only at low T where a similar cancellation can take 
place. Therefore we predict a non-monotonic tempera- 
ture dependence of the net critical current in this case. 

Conclusions. We have found a low temperature 
even/odd sequence of the 0- and tt- states of S/AF/S 
junctions and an additional novel n — transition with 
increasing temperature in odd junctions. The even/odd 
effect is caused by the different symmetries of the even 
versus odd AF interfaces and the corresponding S ma- 
trices, and is revealed in qualitatively different temper- 
ature dependencies of the critical currents for even and 
odd barrier thicknesses. The n — transition with vary- 
ing temperature is induced in the odd junctions by the 
interplay of the spin-split Andreev bound states. 
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